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Abstract. We have computed the charge that develops on an SQN in space as a result of balance 
between the rates of ionization by ambient gammas and capture of ambient electrons. We have 
also computed the times for achieving that equilibrium and binding energy of the least bound SQN 
electrons. We have done this for seven different settings. We sketch the calculations here and give 
their results in the Figure and Table II; details are in the Physical Review D. 79. 023513 (2009). 

INTRODUCTION 

It is now three dozen years since Bodmer's paper [1] suggesting that, while matter 
made of just up and down quarks is not stable against decay into baryons, matter 
made of u, d, and s might be. It is two dozen since Witten's [2] seminal paper 
hypothesizing production of nuggets of strange quark matter (SQM) in the early 
universe when a phase transition might have caused strange quark nuggets (SQNs) to 
condense out of the quark soup. But there has been no experimental or observational 
sign of SQM — in spite of a reasonably vigorous program of searching. It is, however, 
important to pursue as wide a range of tests as possible, both for potential SQN 
discovery and for, lacking that, placing tighter bounds on loci in parameter space 
where SQM could be hiding. 

The basic idea in concluding that quark matter with 3 kinds of quarks is more 
likely to be stable than quark matter with only two is that the potential energy of a 
system with N q quarks is the same independent of what flavor the quarks have. It 
only depends on what colors they have. However the kinetic energy of the system 
depends on both colors and light flavors (i.e., on how many fermi seas). Thus, by the 
Pauli principle, for fixed quark total, the more kinds of flavors, the lower the fermi 
energies and the higher the likelihood of binding. Introducing the c quark doesn't 
help because its mass is greater than the nucleon mass therefore costing more energy 
than is gained by adding another flavor. 

Phenomenologically, the most important feature of SQM is its nuclear density, 
roughly 5 x 10 14 gcm -3 . Also, as Farhi and Jaffe showed [3], if SQM is, in fact, 
stable then, as the number of quarks in the system increases the binding energy per 
quark increases to an asymptotic value. Based on these properties deRujula and 



TABLE I: Some Strange Quark Nugget Searches. References for the data are found in (a)Sandwciss 
[9], (b)Pricc [1], (c) Spiering [11], (d) Banerdt et al. [12], (e) Herrin et al. [13]. 



Experiment/Observation Mass Range (g) Result 

AMS° 1CT 24 - 1(T 22 not done 

RHIC" < 3 x 1(T 21 not found 

Mica Tracks 6 1CT 20 - 1CT 14 << p DM 

ICE CUBE C 1CT 3 - 1(T 2 not done 
Seismometers: 

Future Lunar d 10 3 - 10 6 not done 

Apollo" 10 4 - 10 6 < pdm/W 

USGS Reports c 10 s - 10 s < p DM 



Glashow [4] sketched out what could be learned from a variety of observations and 
experiments if the dark matter (DM) in the galaxy (about 10 _25 gcm -3 ) is SQM in 
the form of SQNs of mass M. Alcock, Farhi and Olinto [5] studied the properties of 
"neutron stars" if they are actually strange quark stars (SQSs). 

A recent development of importance was the recognition that the attractive, fla- 
vor blind quark-quark color force leads to Cooper pairing [6] adding significantly to 
the binding energy of the system. However, on the negative side, experiments and 
observations have lagged far behind theory. They have failed to detect any evidence 
for SQM existence. Some past, present and future SQM searches are summarized 
in Table I. Beside the absence of SQN detection, there are two sets of observations 
that cast doubt on the prediction [5] that neutron stars are actually SQS's. The first 
is the longstanding problem of pulsar glitches. These are sudden changes in pulsar 
frequency, successfully understood on the basis of quakes in the crust of nuclear (not 
quark) matter around the neutron matter core. They are not easily understood on the 
basis of an SQN core, for most models, which can only support a crust too thin for the 
quakes needed [5]; recently, however, theory has managed to produce a "specialized" 
SQN model which is capable of large enough quakes [7], but there does not yet seem 
to be any other confirmable predictions of the model. Another problem with SQS's 
is that from Fourier analysis of quasi periodic oscillations (QPOs) in"superbursts," 
nuclear explosions on soft gamma ray repeaters. The frequencies have been shown [8] 
also to indicate a crust depth greater than the limit found by [5]. 

The object of the present work is to compute the charges that should develop on 
an SQN in space subject to ambient, ionizing radiation and deionizing capture of 
ambient free electrons. As a basic step in devising SQN space searches. 



CALCULATIONS 

We considered seven "settings" with both ionizing radiation producing net charge 
on the SQN and ambient electrons neutralizing that charge. The rate for the first 
falls with increasing (number of electrons ionized; hence total positive charge of 
quark lattice with remaining electrons), while that for the second starts from near 
zero and rises with increasing Z^j. We found the Zn for equality. The seven settings 
are, in the order of decreasing equilibrium charge: (1) the sun shining on an SQN 
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FIG. 1: SQN charge Zn{M). From the top, the 7 curves are in the order described in the text and used for 
the 7 rows of Table II. Mass, M, is in grams. 

near the Earth during a solar X-ray flare; (2) the (extra-galactic) diffuse background 
radiation (DBR) shining on an SQN in intergalactic space; (3) the diffuse background 
radiation shining on and SQN near the Earth; (4) an SQN in the primordial universe 
at recombination; (5) the galactic diffuse radiation (GDR) at the center of the galaxy 
(COG) shining on an SQN located near the COG; (6) the GDR shining on a solar 
system SQN near the Earth; and (7) an SQN in the CMB today (ignoring the DRB 
and other radiation). Our model for the structure of the SQN is the so-called "color- 
flavor locked" model [6] in which there is Cooper pairing of quarks sufficiently strong 
that, in the bulk of the SQN the quarks occur in trios with zero charge, no color and 
antisymmetry in space wave functions. All the net quark charge occurs on the surface 
where there are only u- and d-quarks, which have longer Compton wavelengths than 
s-quarks. The latter must be absent from the SQN surface. This means that, for m s 
in the region of 150 MeV, the charge of the lattice of 3A quarks is about 0.374 2 / 3 .We 
also computed the time to reach equilibrium in the charge (Zn) from r = Z^/Z. The 
results are collected into Table III. 

RESULTS 

Figure 1 shows the results as a function of SQN mass (M). We get these results by 
setting the rate (Z + ) at which the SQN loses electrons due to the radiation equal 
to the rate (ZJ) at which it gains electrons from capture of ambient electrons. This 
yields the equilibrium value of the number liberated (Z^). We did this for SQN 
masses ranging from 10 -20 g to 10 +30 g. The bends in the curves of Fig. 1 occur where 
the Bohr orbit of the most loosely bound electron becomes smaller than the radius of 
the SQN. This, as a function of M, occurs at roughly a nanogram. 

Our expression for the rate of ionization is 

Z + = AntfcF^E' > E) (1) 

where F(E' > E) means all the photons with energy E' greater than E. The rate of 
ionization, Z + , as a function of E, decreases with E (as E increases there are fewer 
photons with greater energy). Binding energy magnitude, Eg, is related to Zn by 



TABLE II: Parameters for the seven curves of Fig. 1. E is in eV, n e has units of cm 3 , v e has units of cm/sec. 
F(E) is for the non-thermal distribution cases. F{E) = AE- 3 - 05 ~//(cm 2 seV) where average temperatures 
are not known, we assume 100 K. 



SQN Location 


Radiation 




Ve/W 5 


Solar X-ray 








Flare at 1 AU 


T = 10 3 eV 
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500 


IGM 


F(E); Log(A)=7.1 


io- 7 


60 


Near Quiet Sun 


T = 0.5 eV 


7 


500 


IGM Pre Recombo 


CMB T = 0.26 eV 


287 


300 


Galaxy Center 


F(E); Log(A) = 11.57 


40 


60 


Near Sun(DRB) 


F(E); Log(A)=7.64 
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Today (CMB) 


T = 2.75K 
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Setting 


M 1/3 T Eq (y) 


E B (eV) 


E fl (eV) 






M > 10- 10 g 


10 g 


Galactic Center 


38 x 10" 4 


4.5 


52 


IGM Today: DBR 


151 


100 


760 


Solar system: 








during X-ray flare 


1.8 x 10" 5 


3.5 x 10 4 


3.9 x 10 4 


from GDR 


1.0 x 10 -5 


0.65 


12 


Quiet Sun 


1.7 x 10" 5 


12 


18 


Recombo with CMB 


2.6 x 10~ 7 


8.0 


10 


Today from CMB 


92 


8.7 x 10~ 3 


0.012 



Eb = Z^a/b where b is the maximum of and cib/Zn with ag the hydrogen Bohr 
radius. We take two kinds of photon spectra: a power law falloff for the diffuse 
radiations of settings 2, 5, and 6 and thermal distributions for the other four. Our 
parameters for the settings are given in Table II. Our expression for electron capture 
is given by 

Z_ = ^r 2 N v e n e {l + / e ) (2) 

where f e is a focusing factor, Safronov effect [14 ] taking into account that, with its 
net positive charge, the SQN attracts electrons so that its cross section for electron 
capture is larger than geometric. f e can be written Eb/E £ with E e simply the ambient 
electron's kinetic energy. Other than this effect, our expressions have only geometric 
cross sections with no reference to the underlying particle physics cross sections. 
This is because the latter, off electrons, are on the order of 10~ 25 while the surface 
density of electrons is closer to 10~ 24 so that the probability is that, if the incident 
electron or photon is within the geometric cross section, the reaction will occur. This 
approximation would need to be revisited for a setting involving relativistic electrons. 
These results are discussed in more detail in Abers et al. [16]. 



DISCUSSION 



Given the results of the previous section, the next question is how can the SQN 
charge in these seven settings be exploited to find the SQNs. We give a few directions 
that might be worth examining. 

1. Direct detection. Direct detection needs a system that can identify a charge to 
mass ratio as predicted by the results in the Figure. We begin with the expected 
frequency, dN ev /dt of events (passages into or through a detector of area A), 

dN eV /dt = n S Q N v SQN A (3) 

where n and v are the number density and speed of the SQNs. Putting in the local 
DM density (0.3 nucleon masses/cm -3 = n divided by M for usqn-) an d the galactic 
virial velocity, gives us that, if we want to instrument an area A and ask how many 
events we should get in a time r, the result is 10~ 17 At / M . If a square kilometer could 
be instrumented, nanogram SQNs could be detected at rates up to 100 each second, 
and one gram ones at about one each year, 

2. Unique line. If an SQN were located close enough to a source of high energy 
gammas whose energy is above twice the electron mass (i.e., over 1.1 MeV), the 
SQN would likely emit gammas of characteristic energy about 1.02 MeV. Madsen 
has studied vacuum breakdown into electron-positron pairs in such a case [15]. The 
magnitude of the binding energy of the least bound electron would be 2m e , if one 
waited long enough, from the balance between energetically favored decay of the 
vacuum into an electron positron pair with the electron being captured and further 
ionization by photons, The time for such vacuum decay increases exponentially as 
the two thirds power of the SQN mass. For a given mass, the time decreases as 
Z^f. Thus there will be a range of masses for which the system will radiate gammas 
of approximately 1.02 MeV. For this range of masses, the SQN will ionize until the 
quark-electron total charge is large enough to make the most loosely bound electron 
bound by 1.02 MeV or more, and large enough to keep the time for vacuum decay of 
the same order as the time for ionization. The 1.02 MeV gammas themselves would 
come from electrons being captured by the SQN at a time when the balance between 
further ionization by ambient gammas and vacuum production of electron-positron 
pairs keeps the binding energy of the least bound electrons at or near that value. 
Calculations of these rates can be done, if needed, with some precision using work of 
three of us solving the Dirac equation for the wave functions with binding energies 
E B = —2m e , see [17]. 

3. SQM in the sun. If SQM is the the lowest energy state of baryonic matter, a fair 
amount of it must pass through the sun, but lose energy during the passage leading 
to capture. Calculating the amount captured needs to take into account the size of 
the charge that develops on the SQN since that charge, through the Safronov effect, 
will enhance the rate of SQN energy loss beyond that of the simple geometric cross 
section usually used in SQN energy loss estimates. Another new feature in computing 
the rate of SQN energy loss is that ionizing the SQN passing through the sun will not 



be governed by electrons being ejected by photons. The sun contains more electrons 
than photons, so the SQN charge will be determined by a balance between ionization 
by fast electrons and deionization by capture of slow electrons. 
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